Determining the carrier-envelope phase of intense few-cycle laser pulses 
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The electromagnetic radiation emitted by an ultra-relativistic accelerated electron is extremely sensitive to 
the precise shape of the field driving the electron. We show that the angular distribution of the photons emitted 
by an electron via multiphoton Compton scattering off an intense (/ > 10'^'' W/cm^), few-cycle laser pulse 
provides a direct way of determining the carrier-envelope phase of the driving laser field. Our calculations take 
into account exactly the laser field, include relativistic and quantum effects and are in principle applicable to 
presently available and future foreseen ultra-strong laser facilities. 
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The rapid development of ultra-fast optics has been al- 
lowing the investigation of physical processes at shorter and 
shorter time scales. Time compression of laser pulses to du- 
rations T below two laser cycles and even down to one cycle 
has been demonstrated in the mid infrared [1] (r « 39 fs at 
a central laser wavelength of A — 12 /im), in the near in- 
frared [2] (t « 4.3 fs at A = 1.55 ^m), in the optical [3] 
(r < 4 fs at A = 0.7 /im) and in the extreme ultraviolet 
(XUV) domain [4] (r w 80 as at A = 12 nm). Laser-matter 
interaction in this "few-cycle" regime shows features, which 
are qualitatively new with respect to the more conventional 
"many-cycle" regime: the response, for example, of atoms 
and molecules becomes sensitive to the precise temporal form 
of the electromagnetic field of the laser and, in particular, to 
its carrier-envelope phase (CEP), i. e. the phase difference 
between the carrier wave and the envelope function (see the 
recent review [5] and the references therein). Viceversa, the 
knowledge and the control of the CEP of a laser pulse allows 
in turn to control physical processes like, for example, atomic 
ionization or above-threshold ionization (ATI). So far, few- 
cycle pulses have been produced with intensities below the 
relativistic threshold, corresponding in the optical domain to 
laser intensities / of the order of 10^* W/cm^. Experimental 
determination of the CEP for few-cycle pulses of intensities 
up to / = lO'^'^-lO^^ W/cm^ has been achieved by a stereo 
ATI measurement technique proposed in [6], with an accuracy 
of about 7r/300 [7]. Other methods employed to measure the 
CEP are attosecond streaking [8] and THz-spectroscopy [9]. 
However, these methods are not applicable for laser pulses of 
intensities above / = 10^^ W/cm^ when relativistic effects 
become increasingly important. 

The generation of intense laser pulses is intimately con- 
nected with temporal compression. Not only because tighter 
temporal compression implies, of course, larger pulse inten- 
sities at a given laser energy and waist size, but also because 
available larger intensities allow potential discovery of new 
physical processes at shorter time scales, opening the possi- 
bility of exploiting them to generate even shorter pulses. In 
the present context a laser pulse characterized by a peak elec- 
tric field £ and by a caiTier angular frequency u, is indicated as 
"intense", if the parameter ^ — \e\£/muj is much larger than 
unity, where e < is the electron's charge, m its mass and 



units with /i = c = 1 are employed, as throughout this work. 
An electron in such an intense laser field is ultra-relativistic, 
i. e. its Lorentz factor is much larger than unity [10]. In 
[11] the possibility of generating single-cycle optical laser 
pulses with peak intensities larger than 10^° W/cm^, corre- 
sponding to ^ > 10, is investigated theoretically. Also, in [12] 
the production of few-cycle, intense XUV bursts is envisaged 
by employing relativistic harmonic generation by a planar tar- 
get. Moreover, the Petawatt Field Synthesizer (PES) laser sys- 
tem under construction in Garching (Germany) aims at opti- 
cal laser intensities of the order of 10^^ W/cm^ (^ w 100) by 
compressing an energy of 5 J to only 5 fs, corresponding to 
less than two laser cycles [13]. Finally, at the Extreme Light 
Infrastructure (ELI) [14] facility, unprecedented laser inten- 
sities of the order of 10^^-10^^ W/cm^ are envisaged, with 
pulse durations of about 10 fs. Therefore it is highly desirable 
to have a procedure to determine the CEP of few-cycle laser 
pulses also when their intensity largely exceeds the relativistic 
threshold. 

In the present Letter we provide a method of determin- 
ing in principle the CEP of an intense (/ > 10^" W/cm^) 
few-cycle laser pulse. The method exploits precisely the spe- 
cific features of the electromagnetic spectrum emitted by an 
ultra-relativistic electron, in particular that it emits radiation 
almost exclusively in a naiTow cone of aperture m/ee ^ 1 
along its instantaneous velocity, with being the electron's 
energy at the emission time [10]. Since intense laser pulses 
are of interest here, the laser field is taken into account exactly. 
Moreover, at such high laser intensities and electron energies, 
quantum effects may play in general a crucial role and they 
are also included by performing the calculations in the frame- 
work of quantum electrodynamics (QED) in the Furry picture 
[15]. This requires the exact solutions of the Dirac equation 
in the presence of the external field as initial and final elec- 
tron quantum states. Therefore, we restrict the quantum case 
to scenarios where the external field can be approximated by 
a plane wave and employ electron VoUcov states [15]. If quan- 
tum effects are negligible our method can be generalized to 
hold in principle for an arbitrary external field and we will 
consider the case of a focused Gaussian beam. We show that 
the angular distribution of the radiation emitted by the elec- 
tron is particularly sensitive to the CEP of the driving pulse 
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and that a theoretical accuracy of about tt/IO can in principle 
be achieved. We also derive an analytical formula connecting 
the angular aperture of the spectrum with the value of the CEP. 

We start by describing the few-cycle laser beam as a pulsed 
plane wave linearly polarized along the x direction and prop- 
agating along the positive z direction. The amplitude and the 
central angular frequency of the pulse are indicated as £ and 
CO, respectively. Then, the electric field £{(j)) of the wave 



depends only on the phase 



uj{t — z), with 



= U!{1, 0, 0, 1) and can be written as 5(0) = Eipgicf) x, 
with tp£ {(j)) being an adimensional function. In order to avoid 
the appearance of unphysical static (dc) components in the 
electric field [16], we choose the electromagnetic vector po- 
tential as «4.((/)) = AipAi'f)) X, where A = —£/co, 

^.w^l'^^'W""^^^^"^ ^"^'''^"^ (1) 

[O elsewhere, 

and = ~u>dA.{(f>)/d(p. Here N is the pulse duration 

in units of the laser period 27r/w, 0o is the CEP of the pulse 
and the sin^-envelope ensures that the electric field and its 
derivative vanish at (f) = and (f> = 2N'k. As we have al- 
ready mentioned, the interaction of the electron and the ex- 
ternal plane-wave field is taken into account exactly in the 
calculations by employing Volkov states as quantum in- and 
out-states of the electron and by working in the Furry pic- 
ture [15]. The interaction between the electron and the quan- 
tized photon field instead scales with the fine-structure con- 
stant aQED = « 1/137 and in the parameter regime 
of our interest here it can be accounted for perturbatively 
up to first order. In this approximation the electromagnetic 
spectrum emitted by the laser-driven electron can be calcu- 
lated from the probability that the electron emits one pho- 
ton. We assume that the incoming electron has spin s and is 
counter-propagating with the laser pulse, therefore its initial 
four-momentum is = (e, 0, 0, —p), with e = y/rnF^\-^. 
The outgoing electron instead has spin s' and four-momentum 
p/fj. _ (g'^p')^ with e' — i/m^ + Finally, the emitted 

photon has four-momentum A;' ^ = (cj', fe'), with a;' = |fe'| 
and its polarization states are described by the four-vectors 
e'Jf, with r' E {1,2}. Then, the transition matrix ele- 
ment Sfi of our process can be cast into the convenient form 
Sf, = (27r)(5 {p'_ + kL' P-) {2tt)^S [p'^ + k'^ - pi_) Mf„ 
where for a general four-momentum = {q'^,C[) the no- 
tation q_ = g° — Qz and q± = {qx,<ly) has been intro- 
duced. The precise expression of the amplitude Mfi is rather 
involved and it is not necessary to report it here. We only 
note that it can be written as M ji = Yl^=o '^ifh where Cj 
are coefficients weakly dependent on the physical parame- 
ters of the problem and where the three functions fj with 
/j. = /^dr?^^(7?)e'/o''dV[aV'.4(V)+/3V'i(V)+7] contain all 
the relevant dynamical information of the process. Here we 
have introduced the important parameters a = —m^k'^/{kp'), 
P = m?^'^k'_/2p_{kp') and 7 = w' (e - pcosi?) /(fcp'), 
with [kp') = kf^p'f^ = oj{e' -p'^) = a;[e + p - a;'(l + cost?)] 
and TT — '& and (p being the spherical angular coordinates of 



the emitted photon, assuming the positive z axis as the polar 
axis. It can be shown that the function /o can be expressed in 
terms of the functions /i and /2 as /o = — (/3/2 + a/i)/!- 
Starting from the above quantity Mfi, one can calculate the 
emitted energy spectrum dS" / dfidw' (average energy emitted 
between co' and oj' + dco', in the soUd angle dO = sin 'ffd'&dip) 
as 



d^ 



dfidw' 167r3 
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It can be seen that, as it must be, the above expression of 
the emitted energy spectrum reduces to its classical counter- 
part when the energy of the emitted photon is much smaller 
than the initial electron energy, i. e. in the Umit uj' <^ e. 
In turn, this occurs if the parameter x = (e + p)£,Lo/m^ 
is much smaller than unity [17]. From a physical point of 
view, the parameter x is the laser's electric field ampUtude 
in units of the QED critical field Ecr = rv? I\e\ in the rest 
frame of the incoming electron. Effects of the laser's pulse 
shape on classical multiphoton Thomson scattering (x <C 
1) at laser intensities around the relativistic threshold have 
been studied in [18-20]. CEP effects have also been investi- 
gated in Schwinger electron-positron pair production in time- 
dependent electric fields [2 1 ] , but, for pair production to occur 
at all in that case, laser intensities have been considered larger 
than 10^7 W/cm2. 

Now, as we have mentioned, we are interested in the ultra- 
relativistic regime in which ^ 1. Also, we will consider 
situations in which ~ e, where the interplay between the 
initial electron energy and the laser intensity produces rich 
dynamics of the electron in the laser field. In this parame- 
ter regime, if x *C 1 (classical case) the electron mainly emits 
frequencies of the order of cj' ^ a;^^, while if x ^ 1 (quantum 
case) the electron mainly emits in the energy range w' ~ e. In 
both cases, one can see that the three parameters a, /? and 
7 appearing in the exponential in the functions /, are all of 
the same order and very large. This implies that the functions 
/j can be evaluated by applying the saddle-point method to 
the integrals in 77. In the classical limit, this circumstance re- 
flects the following physical feature: the spectrum emitted by 
an ultra-relativistic electron along a direction (1?, y>) is mainly 
determined by those parts of the electron trajectory where its 
velocity points along (t?, <^), within a small angle of the or- 
der of m/ce <^ 1, with Ee being the electron's energy at the 
emission time, and the other parts of the trajectory give an ex- 
ponentially small contribution [10]. This is the physical rea- 
son why the energy spectrum emitted by an ultra-relativistic 
electron provides detailed information about the electron tra- 
jectory (in the classical regime) and, in turn, about the precise 
form of the driving external field. As we will show below, this 
last feature remains true also in the quantum regime. In fact, 
from the general theory of the saddle-point method it follows 
in our case that in the region of parameters in which the sad- 
dle points have a large imaginary part, the functions /j are ex- 
ponentially smaller than in those regions in which the saddle 
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points are almost real [17]. If we fix the initial energy of the 
electron, the laser intensity and the energy of the emitted pho- 
ton, the functions fj depend only on the photon emission an- 
gles d and ip. We have shown that the region where the saddle 
points are almost real (their imaginary part scales as 1/^ <C 1) 
corresponds in the classical limit exactly to the angular region 
in 1? and (p, which is spanned by the electron velocity when 
it moves inside the laser pulse. This is fully consistent with 
the mentioned classical circumstance that an ultra-relativistic 
electron emits almost only along its velocity. By consequently 
imposing that the saddle points rjs of the phase in fj are real, 
the following condition must be fulfilled 



< 



P 



tan 



cos ip < 1pA.,i 



(3) 



where V'.A.min {ipA,m-dx) dcnotcs the minimal (maximal) value 
of the function ipAi'P) — Jq ^4'' i'si'P')- The values of V'^.min 
and V'.A.max depend on the CEP 0o and the above equation 
yields the bounding angles for the emission cone as functions 
of the CEP. It is worth observing that, since we work in the 
ultra-relativistic regime e ^ m, then p k, e and Eq. (3) pro- 
vides universal conditions depending (apart from the shape of 
the laser field, of course) only on the ratio e/m^. 

The above considerations equally apply in the classical and 
in the quantum regime and the spectrum (2) in the limit x *C 1 
goes to the classical expression derived by first solving the 
Lorentz equation and then plugging the resulting trajectory 
into the Lienard-Wiechert potentials (see Eq. (66.9) in [10]). 
However, the classical formulas hold in principle for an arbi- 
trary external field. Moreover, it can be shown that the above 
equation (3) can also be employed for an external potential of 
the form ^(0; x, y, z) = AiIja{'1>\ x, y, z) x slowly-varying 
with respect to x, y and z, by simply substituting the function 
'4ja{4') with il;j{{(f>; x,y, z). In Fig. 1 we show the depen- 
dence of the angular emission region in as a function of 
the CEP, obtained from Eq. (3) in the two cases — 2e 
(Fig. la) and — e/7.5 (Fig. lb). In both cases N — 2. 
In view of the first numerical example worked out below, in 
Fig. la we consider a Gaussian beam with carrier wavelength 
A = 1.2 jim (u — 1 eV), spot radius w = 2 fim at zero order 
in the ratio w/zj, < 1, with Zj. — uiw'^ /2 being the laser's 
Rayleigh length [23] (we will see that our formula (3) with 
■0^(0) ~^ i'A{4'\ ^1 U: z) works well also at such tight focus- 
ing). In Fig. lb the case of an external plane wave is consid- 
ered. The radiation is confined around the azimuthal angles 
tf — Q and ip — tt within a small aperture angle of roughly 
A(p ~ m/e <C 1, therefore we consider only the directions 
exactly at = and at ip — tt, corresponding to d > and 
?9 < in Fig. 1, respectively. The figure shows a one-to-one 
dependence of the emission range A-d = ?9max ^ ^^min on the 
CEP (po i'dniin and dm^^ are the bounding angles obtained from 
Eq. (3)). From Eq. (3) one also sees that the most convenient 
range of parameters in terms of accuracy in the determina- 
tion of the CEP is at ^ e (if m£, ^ e, then the electron 
emits almost exclusively into a cone along the laser propa- 
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FIG. 1 : Analytical angular emission range as a function of the CEP 
obtained from Eq. (3) for a two-cycle pulse with = 2e (part a)) 
and — e/7.5 (part b)). In part a) the external field is a Gaussian 
beam with carrier angular frequency cj = 1 eV and spot radius w — 
2 fim. In part b) the external field is a plane wave. 



gation direction with an angular aperture of the order of 1/^ 
independently of the CEP). 

In order to show quantitatively the features of our method 
we consider below two examples. In the first example, we aim 
to investigate a rather realistic situation in which the external 
field is modeled as a focused Gaussian beam also including 
longitudinal field components up to first order in w/zr [23]. 
By recalling the laser's parameters envisaged at the PES [13], 
we set w = 1 eV, N = 2 (corresponding to r = 8 fs), 
w = 2 ^m and ^ = 100 (/ = 10^^ W/cm^). Moreover, it can 
be seen that for the electron densities of beams obtained via 
laser wakefield acceleration [24], the relevant (high-energy) 
part of the spectrum results essentially from incoherent emis- 
sion. We consider an electron beam with a three-dimensional 
Gaussian spatial distribution with waists We = We.y = 5 /im 
and We,z ~ 8 ^m and with a Gaussian energy distribution 
with central energy e = 26 MeV (such that « 2e like 
in Fig. la) and waist We,€ such that We.e/e = 2 % [24]. 
Therefore, since x 2 x 10~^, we can calculate the spec- 
trum by employing the classical formula valid for a single 
electron [10] and then averaging it over the electron distri- 
bution (in the numerical example we have considered a beam 
with Ne — 300 electrons and ensured that our results are not 
significantly altered by increasing N^)- In this example (see 
Figs. 2a and 2b) we aim to show the sensitivity of our method 
and we show two energy spectra for the two different CEPs 
<j>o = — tt/IO (Fig. 2a) and 0o = — 7r/5 (Fig. 2b). The 
white horizontal lines indicate the values of "i^min and 'dmax as 
predicted by Eq. (3) generalized to the case of a Gaussian 
beam. The plots show the very good agreement between the 
analytical predictions and the numerical results. Moreover, 
it is evident that the angular aperture of the emission region 
is very sensitive to the CEP: a small change |A0o| — tt/IO 
in the CEP changes the minimum (maximum) emission an- 
gle by about 3° (5°). In Figs. 2c and 2d we consider another 
example, in which an ultra-strong attosecond XUV pulse, as 
those theoretically envisaged in [12] is employed (oj — 50 eV, 
T = 160 as and ^ = 20, which can be obtained if the field is 
focused to about w = 100 nm) and an electron with initial 
energy of 75 MeV (corresponding to 7.5 as in Fig. lb). In 
this case, x = 0.6 and quantum effects are already important. 
Since here A = 25 nm and w = 100 nm, the quantum results 
valid in the plane-wave approximation should apply with suf- 
ficient accuracy. With these parameters the method turns to be 
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FIG. 2: (Color online) Energy emission spectra dS' /dflduj' in sr^^ 
via Eq. (66.9) in [10] (parts a) and b)) and via Eq. (2) (parts c) and 
d)) for the two sets of parameters described in the text. In parts a) 
and b) quantum effects are negligible and it is (fio = — tt/IO (part 
a)) and ijio = — 7r/5 (part b)). In part c) and d) quantum effects are 
important and it is <^o =0 (part c)) and <^o = 7i"/4 (part d)) (the 
almost vertical red line indicates here the quantum cut-off frequency 
uj'j^j = (e + p)/(l + cos ■£))). The horizontal white lines indicate the 
boundary of the emission range determined analytically from Eq. (3) 
generalized to the case of a Gaussian beam for parts a) and b). 



less precise as compared to the above example and a change in 
the CEP of about tt /4 produces a change in the angular aper- 
ture of a few degrees (in Fig. 2c it is (/)o = while in Fig. 2d 
4'o = '^/^)- In this case we observe an excellent agreement 
between the numerical and the analytical values of i^min and 
i^m-dx- Moreover, the almost vertical red line shows the posi- 
tion of the cut-off emission frequency cj^j determined from 
the analytical formula = {e + p)/{l + cos'd) [17]. This 
is a typical quantum effect: due to energy-momentum conser- 
vation, the electron cannot emit a photon at an angle ■& with 
energy larger than or equal to o;^^. Therefore, also the quan- 
tum cut-off frequency uj'j^j is affected by the CEP, through the 
emission aperture Ai?. Finally, we note that the bright features 
of the spectra in Figs. 2c and 2d along the direction d = can 
be explained classically: in a relatively large part of the tra- 
jectory inside the laser field the electron velocity is observed 
to point just along that direction. 

Experimental uncertainties in the laser intensity may alter 
our results. The intensity of a strong optical laser beam can be 
measured nowadays with a relative uncertainty A/ //of about 
10 % [22]. If one includes a coiTesponding uncertainty A^ in 
the value of ^ in Eq. (3), one obtains that the induced uncer- 



tainty A?9min/max in the predicted values i?min/max is approx- 
imately given by AlJmin/max « 4rne|V'.A,inin/max|A^/(4e^ + 

™^C^I'0^,min/maxP)- In the example in Figs. 2a and 2b we ob- 
tain At?n,in/[nax ~ 2.8°. In the second example, it is difficult to 
estimate the value of A/ //. However, since e w 7.5 then 
Ai9,Tiin/max ~ ^nd even an uncertainty of about 

50 % in the intensity is acceptable. Thus in both cases we can 
conclude that these uncertainties do not conceal the effect of 
the CEP. By repeating the simulation with a different temporal 
pulse shape (a Gaussian one) we have observed alterations to 
the spectra, much smaller than those due to the uncertainty in 
the laser intensity. We finally note that for pulses comprising 
more than three laser cycles the discussed effect is too small 
to be significant. On the other hand, the CEP effect is even 
larger than here for pulses including only one laser cycle. 
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